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Abstract: Recent molecular cytogenetic studies demonstrate that extensive centromere-telomere fusions are the main 
chromosomal rearrangements underlying the karyotypic evolution of extant muntjacs. Although the molecular mechanism 
of tandem fusions remains unknown, satellite DNA is believed to have facilitated chromosome fusions by non-allelic 
homologous recombination. Previous studies detected non-random hybridization signals of cloned satellite DNA at the 
postulated fusion sites on the chromosomes in Indian and Chinese muntjacs. But the genomic distribution and 
organization of satellite DNAs in other muntjacs have not been investigated. In this study, we have isolated four satellite 
DNA clones (BMC5, BM700, BM1.1k and FM700) from the black muntjac (Muntiacus crinifrons) and Fea’s muntjac (M. 
feae), and hybridized these four clones onto chromosomes of four muntjac species (M. reevesi, M. crinifrons, M. 
gongshanenisis and M. feae). Besides the predominant centromeric signals, non-random interstitial hybridization signals 
from satellite I and II DNA clones (BMC5, BM700 and FM700) were also observed on the arms of chromosomes of these 
four muntjacs. Our results provide additional support for the notion that the karyotypes of M. crinifrons, M. feae and M. 
gongshanensis have evolved from a 2n = 70 ancestral karyotype by a series of chromosome fusions. 
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Muntiacus have experienced rapid апа radical 
chromosome evolution within the last two million years 
(Fontana & Rubini, 1990; Yang et al, 1995; Wang & Lan, 
2000). The extreme diversification in chromosomal 
number and structure make them an ideal model for 
studying karyotypic evolution. Among the five extant 
muntjacs that have karyotype reports, the Indian muntjac 
(M. muntjak vaginalis) has the lowest chromosome 
number known in mammals, 2n = 6 ($) and 2n = 7 (2) 
(Wurster & Benirschke, 1970); the Chinese muntjac (M. 
reevesi), a close relative of the Indian muntjac, has a 
karyotype of 2n — 46, all chromosomes are acrocentric 
(Wurster & Benirschke, 1967); the other species have 
intermediate karyotypes: 2n = 8 (©) and 2n = 9 (5) in 
black muntjac (M. crinifrons) (Shi, 1983) and Gongshan 
muntjac (M. gongshanensis) (Shi & Ma, 1988); 2n = 
12,13,14 ($) and 2n = 14 (£) in Fea's muntjac (M. feae) 
(Soma et al, 1983, 1987; Tanomtong et al, 2005). 

Several hypotheses were proposed to explain the 
great karyotypic diversity in Muntiacus. The most 
well-known hypothesis is the tandem fusion hypothesis, 
first proposed by Hsu et al (1975), suggesting that the 2n 
= 69/72 karyotype of M. muntjak vaginalis could have 
evolved from 2n = 46 M. reevesi-like ancestral karyotype 
through extensive tandem fusions and several centric 
fusions. Subsequent studies, including conventional 
comparative cytogenetics (Shi et al, 1980; Elder & Hsu, 
1988; Fontana & Rubini, 1990), chromosome painting 
(Yang et al, 1995, 1997a; Yang, 1998; Chi et al, 2005a; 
Huang et al, 2006a), FISH mapping of centromeric 
sequences, telomeric sequences, cosmid clones and BAC 
clones (Scherthan,1990; Lin et al, 1991; Lee et al, 1993; 
Froenicke et al, 1997; Li et al, 2000a; Hartmann & 
Scherthan, 2004; Chi et al, 2005b; Huang et al, 2006b), 
and combined chromosome painting and satellite DNA 
mapping (Scherthan, 1995; Yang et al, 1997b,d; 
Froenicke & Scherthan, 1997; Huang et al, 2006c), 
provided direct evidence for the tandem fusion 
hypothesis. The chromosomal mechanism underlying the 
karyotype evolution in Muntiacus is well-established 
now: the karyotypes of all extant muntjacs have evolved 
from a common ancestor with a 2n — 70 acrocentric 
karyotype by extensive centromere-telomere fusions and 


several centric fusions. 
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Nevertheless, the molecular mechanism that 
triggered such extensive tandem chromosomal fusions 
remains unclear. Some studies suggested that the 
repetitive DNA families at or near the centromeric and 
telomeric regions might facilitate illegitimate 
recombination between non-homologous chromosomes 
of muntjacs (Brinkley et al, 1984; Bogenberger et al, 
1985, 1987; Benedum et al, 1986; Lin et al, 1991, 2004; 
Lee et al, 1994, 1997; Scherthan, 1995; Lee & Lin, 1996; 
Yang et al, 1997b; Li et al, 2000a, b, 2002; Hartmann & 
Scherthan, 2004). At present, four satellite DNA families 
are found in Muntiacus: satellite DNA families I, П, IV 
and V (Bogenberger et al, 1985; Lin et al, 1991, 2004; Li 
et al, 2000b, 2005). FISH mapping demonstrated that 
satellite I DNA of M. muntjak vaginalis, M. reevesi and 
M. reevesi micrurus and satellite II DNA of M. muntjak 
vaginalis are localized at both the centric regions and at 
non-random interstitial sites along the arms of the 
"fusion" chromosome (Lin et al, 1991, 2004; Li et al, 
2000b). The findings of telomeric repetitive sequences 
present at several interstitial locations іп M. muntjak 
vaginalis chromosomes (Lee et al, 1993; Scherthan, 1995) 
as well as the comparative mapping of the satellite I and 
II families (Li et al, 2000b) indicated that during tandem 
fusions the chromosomal breakpoints localized at 
satellite II DNA regions and subtelomeric regions of the 
ancestral chromosomes. 

Until now, all satellite DNA families of the genus 
Muntiacus were isolated from the genomes of M. 
muntjak vaginalis, M. reevesi and M. reevesi micrurus (a 
subspecies of M. reevesi, 2n — 46), as did the FISH 
mapping of such satellite DNA. To date, no satellite 
DNA has been isolated from the genomes of M. 
crinifrons, M. gongshanensis and M. feae. Here we have 
cloned and characterized four satellite DNA clones from 
the genomes of M. crinifrons and M. feae, and studied 
the chromosome distribution of these four satellite DNAs 
in M. crinifrons, M. feae, M. gongshanensis and M. 
reevesi. In doing so, we hope to provide additional 
insights into the molecular mechanism of the tandem 
fusions that led to the formation of the karyotypes of M. 


crinifrons, M. gongshanensis and M. feae. 


1 Materials and Methods 
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1.1 Cell culture and chromosome preparation 

Fibroblast cell lines of a male M. crinifrons 
(KCB200004), a male M. feae (KCB 91006), a male M. 
reevesi (KCB91001) and a female M. gongshanensis 
(KCB 88003) were obtained from the Kunming Cell 
Bank (Kunming, Yunnan, PR China). Cells were cultured 
at 37'C under 5% CO; in DMEM medium supplemented 
with 10% fetal bovine serum, 100 U/mL penicillin and 
100 релт, streptomycin. Chromosome preparations 
were prepared as previously described (Yang et al, 1995). 
1.2 Molecular cloning of repetitive DNA sequences 

Genomic DNA of male M. crinifrons and M. feae 
were extracted from the liver tissue and cultured 
fibroblast cells of these two species by using a Genomic 
DNA Extraction kit (Bioteke, China). Genomic DNAs of 
M. crinifrons and M. feae were digested with five 
restriction endonucleases: EcoRI, Веі, BamHI, HindIII 
and Xbal, respectively. The digestion products were 
subsequently fractionated on a 1.2% agarose gel. The 
predominant DNA fragments with a size of about 400bp 
were purified and cloned into pUC19 vector (Takara). 
PCR was carried out in a 50uL reaction volume using 
100ng template DNA, Ex-Taq DNA polymerase (Takara), 
and 0.2 umol/L of Satellite II primers (SatII-fw[5'-GAG- 
СТСССТСАСАСАСТСС-37 and SatII-rv [5'-CAG- 
AGCCGACCTAGGATCAC-3']) as previously described 
(Li et al, 2000b). PCR products were fractionated on 
1.596 agarose gel and predominant DNA fragments were 
purified, and cloned into pMD18-T Vector (Takara). 
Clones containing repetitive DNA inserts were screened 
by PCR with the M13 forward and reverse primers. 
1.3 Nucleotide sequencing and analysis 

The satellite DNA clones were sequenced and 
deposited in NCBI GenBank. The accession numbers for 
ВМС5, BM700, ВМІЛК and FM700 аге EU644506, 
EU644507, EU644508 and EU644509 respectively. 
BLASTN 


performed to find similarity sequences in the Database. 


(http://www.ncbi.nlm.nih.gov/blast ^ was 


Tandem Repeats Finder Program (Benson, 1999) was 
used to reveal tandem repeats in the repetitive DNA 
fragment. Dot Matrix and multiple-alignment were 
performed by using DNAMAN software (Version 4). 
1.4 Fluorescence iz situ hybridization (FISH) 

DNA fragments of the satellite DNA clones were 
labeled by nick translation with  biotin-14-dCTP 
(Invitrogen) or FITC-12-dUTP (Roche) FISH was 


performed as described previously (Yang et al, 1997b). 
Biotin-labeled probes were visualized using Cy3-avidin 
(1:1000, Amersham), while FITC-labeled probes were 
detected with a layer of rabbit-anti-FITC IgG (1:200, 
Invitrogen) followed by a layer of FITC conjugated 
goat-anti-rabbit IgG (1:250, Vector Laboratories). After 
detection, slides were mounted in Vectashield medium 
with DAPI 


Laboratories) and covered with 22 mm x 32 mm 


(4'6-diamidino-2-phenylindole, Vector 
cover-slips. FISH images were captured using the Genus 
system (Applied Imaging Corp.) as previously described 
(Yang et al, 2000). Hybridization signals were assigned 
to specific chromosomes or chromosomal regions as 
defined by enhanced DAPI-banding patterns (Yang et al 
1995) or pre-captured G-banding patterns. 


2 Results 


2.1 Molecular cloning and characterization of sate- 
lite I, П, IV DNA clones from the genomes of 
Muntiacus crinifrons and M. feae 
One clone, ВМС5, was isolated from BamHI digests 

of M. crinifrons genomic DNA. Three repetitive DNA 
clones were generated from PCR amplifications of M. 
crinifrons and M. feae genomic DNA using primer 
sequences derived from a white tailed deer satellite П 
DNA sequence. Among them, two clones were isolated 
from M. crinifrons, BM700 and BMI.1k, and one clone 
was isolated from M. feae, FM700. These four clones 
were characterized by sequencing and aligning with the 
known cervid satellite DNA, and shown to belong to 
cervid satellite DNA clones. 

BMCS had a length of 437 bp and 51.49% GC 
content. A BamHI and EcoRI digested DNA fragment, 
with a length of 221bp from nucleotides 1—222 of BMC5 
clone, had 93% homology with C5 clone (M. reevesi 
satellite I clone) with the nucleotides 561—783 (Lin et al, 
1991), and nucleotides 235-437 of ВМС5 clone had 
89% similarity with the nucleotides 20—221 of C5 clone 
(Fig. la). When BMCS5 clone was aligned with FM-sat I 
clone (M. reevesi micrurus satellite I clone) (Lin et al, 
2004), 88% similarity was found between nucleotides 
1-437 of BMC5 clone and nucleotides 775-1211 of 
FM-sat I clone, and there was also a 91% similarity 
between nucleotides 26-437 of BMCS clone and 
nucleotides 4—416 of FM-sat I clone (Fig. 1b). These 
results suggest that BMCS belongs to the satellite I DNA 
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Comparison of satellite I DNA sequences among muntjac species 


a: Schematic diagram of DNA sequences comparison between Muntiacus crinifrons satellite I clone (ВМС5) and Chinese muntjac 


satellite I clone (C5) (EcoRI sites indicated by E and BamHI sites by B); b: Schematic representation of DNA sequences comparison 


between M. crinifrons satellite I clone (BMC5) and Formosan muntjac satellite I clone (FM-sat I). 


BM700 and FM700 are 700-bp PCR products 
amplified from M. crinifrons and M. feae genomes 
respectively using satellite II primers (Li et al, 2000b) 
(Fig. 2a). BM700 is 667bp in length and had 62.37% GC 
content, and FM700 had a length of 660bp and 63.03% 
GC content. Nucleotides 242—634 of BM700 clone had 
84% similarity with nucleotides 141—533 of FM-sat II 
clone (M. reevesi micrurus satellite II clone) (Lin et al, 
2004), and nucleotides 1-195 had 74% homology with 
nucleotides 619-817 of FM-sat IL, and nucleotides 
242-543 of BM700 was of 81% similarity with 
nucleotides 843-1143 of FM-sat II clone (Fig. 2b). When 
FM700 clone was aligned with FM-sat II clone, 81% 
similarity was found between nucleotides 83—634 of 
FM700 and nucleotides 1-555 of FM-sat II clone, and 
nucleotides 1-520 of FM700 had 80% homology with 
621-142 of FM-sat II (Fig. 2c). 
Multiple-alignment was performed among ВМ700, 
FM700 and MMV-0.7 (M. muntjak vaginalis satellite II 
clone) (Li et al 2000b), and 85.69% homology was found 
among these clones (Fig. 2d). These findings indicate 
BM700 and FM700 belong to the satellite II DNA 
family. 

ВМІЛЕ was a PCR product amplified from М. 


nucleotides 


crinifrons genome using satellite II primers (Li et al, 
2002) (Fig. 2a). BM1.1k had a length of 1.1-kbp with 
44.27% GC content. Multiple-alignment was performed 
among BMI.Ik, MMV-1.0 (M. muntjak vaginalis 
satellite IV), MR-1.0 (M. reevesi satellite IV) (Li et al, 
2002), and FM-sat IV (M. reevesi micrurus satellite IV) 
(Lin et al, 2004), 98.17% homology was found among 
these clones (Fig. 2e), suggesting that BM1.1k belongs to 
the satellite IV DNA family. No tandem repeats were 
found in these repetitive DNA elements using the 
Tandem Repeats Finder Program (Benson, 1999) and no 
internal sub-repeats were detected by dot matrix analysis. 
2.2 Chromosome distribution of satellite I (BM C5), 

satellite II (BM700) and satellite IV (BMI.1k) 

DNAs in M. crinifrons 

The probe from M. crinifrons satellite I element 
(BMC5) was hybridized to the metaphase chromosomes 
of a male M. crinifrons (2n = 92). Predominant signals 
were observed at the centric regions of all chromosomes 
except for the Y chromosome which is relatively weaker, 
but non-random signals were also found at some 
interstitial locations along the chromosomes of M. 


crinifrons, with seven interstitial signals on Chrl, 
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Fig. 2 Molecular cloning and characterization of sate-llite I, П, IV DNA clones from the genomes of 


a: Electrophoretic analysis of PCR products. PCR products were amplified with a pair of Satellite II primers (Li et al, 2000b, 2002) from 
three muntjac species (IM, Muntiacus muntjak vaginalis; BM, M. crinifrons; Feas, M. feae) and ‘con’ stands for negative control without 
template. The products were fractionated on a 1.5% agarose gel. Three bands of 0.7, 1.1 and 1.4kb are detected in these species; b: 
Schematic illustration of sequence comparison between M. crinifrons satellite II DNA clone (BM700) and Formosan muntjac satellite II 
clone (FM-satII); c: Schematic diagram DNA sequences comparison between M. feae satellite DNA II clone (FM700) and Formosan 
muntjac satellite П clone (FM-satII); d: Multiple-alignment of M. crinifrons satellite II clone (BM700), M. feae satellite II clone (FM700) 
and M. muntjak vaginalis satellite II clone (MMV-0.7) shows 85.69% homology (identical nucleotide sequences from these clones are 
shown in black boxes); e: Multiple-alignment of M. crinifrons satellite IV clone (BM1.1k), MMV-1.0 (M. muntjak vaginalis satellite IV) 
(Li et al, 2002), MR-1.0 (M. reevesi satellite IV) (Li et al, 2002), and FM-sat IV (M. reevesi micrurus satellite IV) (Lin et al, 2004) shows 


98.1796 similarity. 
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Fig. 3 FISH mapping examples 


a: Satellite I DNA (ВМС5) signals (red) are detected at centric regions of all Muntiacus crinifrons chromosomes except Y chromosome 
as well as at non-random interstitial fusion sites; b: Cy3 image of the same metaphase; c: Hybridization of M. crinifrons satellite II DNA 
clone probe (BM700) on a male M. crinifrons metaphase. Arrows indicate the satellite II DNA signals; d: Co-hybridization of satellite II 
(BM700, green) and satellite IV (BM1.1k, red) probes on a male M. crinifrons metaphase; e: Localization of M. crinifrons satellite I 
clone (BMC5) probe on male M. reevesi metaphase; f: Localization of M. crinifrons satellite I (ВМС5, green) and satellite II (BM700, 
red) DNA probes on male M. reevesi metaphase. Satellite II DNA is located distal to satellite I DNA (indicated by arrows). 
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ChrX-4 and Chr1P+4, six interstitial signals on Chrlq 

and Chr2, five interstitial signals on Chr3, respectively 

(Fig. 3a, b). The hybridization patterns were summarized 

on a M. crinifrons idiogram (Fig. 4a). In total, 49 

site-specific autosomal interstitial signals were detected 

by BMCS probe on the nine M. crinifrons chromosomes. 

The probes from satellite II element (BM700) and 
satellite IV element (BM1.1k) of M. crinifrons were also 
hybridized onto the metaphase chromosomes of M. 
crinifrons. Besides the predominant pericentromeric 
signals, several non-random interstitial signals were 
found along the chromosomes of M. crinifrons by 
BM700 probe (Fig. 3c). The probe of BMI.Ik gave 
signals exclusively at Ше centromeres of all 
chromosomes in М. crinifrons. Two-color FISH, 
simultaneously hybridizing satellite II DNA (BM700) 
probe and satellite IV DNA (BMI.1k) probe onto the 
metaphases of male M. crinifrons, demonstrated that the 
hybridization signals of satellite IV (red) co-localized 
with the satellite II signals (green) at the centric regions 
of all M. crinifrons chromosomes (Fig. 3d). 

2.3 Chromosome distribution of satellite I (BMC5) 
and satellite II (BM700 and FM700) DNA in M. 
reevesi, M. feae and M. gongshanensis 
The probe from M. crinifrons satellite I element 

(BMC5) was also hybridized to the metaphase 

chromosomes of M. reevesi (2n = 46), M. feae (2n = 14) 

and M. gongshanensis (2n = 8). Besides the centromeres, 

ВМС5 probe also gave signals at specific interstitial 

locations along chromosome 1—5, 11 of M. reevesi (Fig. 

Зе). The hybridization patterns of BMC5 on the 

chromosomes of M. feae and M. gongshanensis were 

very similar to that of BMC5 on the chromosomes of M. 

crinifrons. Some non-random interstitial signals were 

observed along the chromosomes of M. feae and M. 

gongshanensis. Altogether 28 and 46 site-specific 

autosomal interstitial signals were detected in the 
genomes of male M. feae and female M. gongshanensis 

(Fig. 4b, c). 

To confirm the location of satellite I DNA and 
satellite II DNA on chromosomes of M. 
two-color FISH was performed by simultaneously 
hybridizing the satellite I (BMC5) probe and satellite II 
(BM700 or FM700) probe onto chromosomes of M. 
reevesi. The results indicate satellite II element (BM700, 
red) is distal to satellite I (BMCS, green), and interstitial 


reevesi, 


(ВМС5, also observed оп 


chromosomes 1—5 and 11 of M. reevesi (Fig. 3f). Similar 


signals green) were 
chromosome distribution patterns were also detected on 
the chromosomes of M. reevesi using M. crinifrons 
satellite I DNA (BMC5) and M. feae satellite П DNA 
(FM700) probes (data not shown). 


3 Discussion 


We have successfully cloned and characterized four 
centromeric satellite DNA clones, which belong to three 
different satellite DNA families: I (BMC5), II (BM700 
апа FM700), and IV (BM1.1k), from M. crinifrons and 
M. feae. Mapping these four satellite DNA clones onto 
the chromosomes of M. reevesi, M. crinifrons, M. feae 
and M. gongshanensis allowed us to further investigate 
the distribution and organization of satellite DNAs in 
other muntjac species and the nature of the chromosome 
fusions that lead to the origin of diverse karyotypes in 
different muntjac species. 

3.1 Satellite I clone-BMC5 

The hybridization pattern of BMC5 to the 
chromosomes of M. reevesi closely resembles that of the 
C5 clone probe in the metaphase of M. reevesi (Yang et 
al, 19974; Li et al, 2000b), indicating that BMCS5 satellite 
I DNA, like the C5 centromeric satellite DNA, was 
probably inherited from the ancestral acrocentric 
chromosomes, and represents the remnant of centromeric 
heterochromatin of ancestral chromosomes after tandem 
chromosomal fusion. 

Recent chromosome painting and BAC mapping 
studies demonstrated that 27 and 28 centromere-telomere 
tandem fusions are needed to “reconstruct” the haploid 
karyotypes of M. feae (2n = 14), M. crinifrons and M. 
gongshanensis (both 2n = 8,9) respectively from a 2n = 
70 ancestral karyotype (Yang et al, 1997c; Huang et al, 
2006c). ВМС5 clone probe revealed 49, 28 and 46 
interstitial hybridization signals in the diploid cells of M. 
crinifrons, M. feae and M. gongshanensis, respectively 
(Fig. 4). Moreover, these interstitial hybridization signals, 
apparently mapped to the putative tandem fusion sites 
along the chromosomes of M. crinifrons, M. feae and M. 
gongshanensis defined previously by chromosome 
painting and comparative BAC mapping. Although 28 
interstitial hybridization signals detected by ВМС5 clone 
probe in M. feae was only half of the number of putative 
fusion sites (54), 49 and 46 interstitial hybridization 
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Fig.4 Summary of the non-random interstitial signals of Muntiacus crinifrons satellite I clone (BMC5) on a high- 


resolution G-banded idiogram of M. crinifrons (a), M. feae (b) and M. gongshanensis (c) 


The ideograms were modified from Yang et al, 1998. The homologous chromosomes or segments of the Chinese muntjac (MRE) 


indicated on the right of M. crinifrons, M. feae and M. gongshanensis chromosomes. The interstitial sites are indicated by *. 
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signals in M. crinifrons and M. gongshanensis detected 
by BMCS clone probe respectively were very close to the 
number of putative fusion sites (56). The absence of 
interstitial hybridization signals in some putative fusion 
sites of M. feae 1s most likely due to variations in copy 
number of satellite I DNA (if the copy number is too low, 
the hybridization signals will be too weak to be detected 
by FISH). These data, together with the interstitial 
signals detected on M. reevesi chromosomes 1—5 and 11 
by ВМС5 clone probe, 


evidence for the origin of M. crinifrons, M. feae and M. 


provide direct molecular 
gongshanensis karyotypes from an ancestral karyotype 
(2n = 70). 
3.2 Satellite II clones - BM700 and FM700 

The FISH mapping results of M. crinifrons (BM700) 
and M. feae (FM 700) demonstrated that satellite II DNA 
was mainly found in the centromeric regions of the 
chromosomes of M. crinifrons (Fig. 3c) and M. feae (data 
not shown). Several interstitial hybridization signals 
were also detected on the arms of chromosomes of M. 
crinifrons (Fig. 3c) and M. feae (data not shown) by 
these two satellite II DNA clone probes. But the number 
of interstitial hybridization signals detected by satellite II 
DNA clone probes was far less than that by satellite I 
DNA clone probes. In contrast, Li et al (2000b) showed 
that the M. muntjak vaginalis satellite II clone (MMV-0.7) 
gave signals on most sites of tandem fusions along M. 
muntjak vaginalis chromosomes. This discrepancy could 
be due to either the variation in the amounts of target 
satellite II DNA between species, or most likely, the 
variation in FISH protocols as the MMV-0.7 probes 
made by us using the same PCR primer reported by Li et 
al (2000b) only gave several interstitial hybridization 
signals in the Indian muntjac genome (data not shown). 
The two-color FISH results of satellite I (BMC5) and II 
(BM700 and FM700) DNA probes further indicate that 
satellite II DNA localized distal to satellite I DNA at 


centromeric regions in M. reevesi (Fig. 3f). This finding 
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